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^ A complete and rigorous determination of the possible ground states for D-wave pairing Bose 

condensates is presented, using a geometrical invariant theory approach to the problem. The order 
parameter is argued to be a vector, transforming according to a ten dimensional real representation 
of the group G —O^iSiVi x (T). We determine the equalities and inequalities defining the orbit 
space of this linear group and its symmetry strata, which are in a one-to-one correspondence with 
the possible distinct phases of the system. We find 15 allowed phases (besides the unbroken one) , 
P3 . with different symmetries, that we thoroughly determine. The group-subgroup relations between 

O ' bordering phases are pointed out. The perturbative sixth degree corrections to the minimum of a 

^ , fourth degree polynomial G-invariant free energy, calculated by Mermin, are also determined. 

U 

Oh' PACS 74.20.De, 02.20.-a, 64.60.-i, MSC 14L24 

'j^ ; I. INTRODUCTION 

Superfluidity and superconductivity are justified on the basis of the macroscopic condensation of Bose quasi- 
' ^ particles. The classical Bardeen, Cooper and SchriefFer theory for superconductivity dates 1957. Soon after, a 
BCS-type transition was proposed for the Fermi system '^He by Anderson and Morel ||l|. Cooper pair formation was 
thought to occur in an L ^ state, to take into account the hard core nature of "^He atom interaction. The superfluid 
phases were actually observed ||], and the nature of p-wave pairing is now well established for '^He. The theory of 
L y^ 0-superfluids is relevant for "unconventional" superconductivity too. The high temperature superconducting 
^►^ (HTS) oxides are anomalous in their non- Fermi liquid normal state properties and in many cases share with heavy 

^— V , fermion superconductors unconventional d-wave pairing behaviour MM- In fact, experiments probing the phase and 
QQ ' the nodes of the gap show a sign reversal of the order parameter, compatible with the d-wave scenario. Nevertheless, 
f^ , other experiments on the same HTS compounds point in the direction of a significant s-wave component. In the 
0^ ' presence of external magnetic fields, defects and interface phenomena there are also evidences of a mixed pairing 
^D I symmetry (see e.g., the Introduction of |g] and references therein). Such a controversial experimental situation is 
^--^ . reflected in the fact that the underlying microscopic mechanism inducing superconductivity in HTS materials is still 
unclear and is one of today's major challenges ^,0]. 

Such a situation has motivated the efforts at studying the macroscopic properties of unconventional superconductors 
through the Landau theory of phase transitions [|| and the identification of the order parameter symmetry has become 
a preliminary task in the construction of viable models describing the attractive nature of the pairing interaction. 
^-^ • Direct experimental evidence of an order parameter unconventional structure lies on multiple phase transitions. The 

Ch \ heavy fermion compounds \J i^x'^^x^^is display four different superconducting phases in the T-x phase diagram and 
O ■ UPta displays three superconducting phases in a T-H phase diagram |^. Furthermore, from power-law temperature 
behaviour of thermodynamic and transport properties (e.g. specific heat, magnetic penetration depth), a non-trivial 
node structure for the gap function may be inferred, compatible with high i-pairing. That is the case also for HTS 
oxides, for which, at present, clear proofs are lacking of the existence of more than one superconducting state. 
, . , Actually, it must be pointed out that, unlike superfiuid '^He, which is an isotropic fermion system, Bloch electrons 

^ • in a superconductor crystal lattice exhibit, in general, a reduced finite symmetry; in addition, the existence of 
imperfections, even in the cleanest samples, can partly destroy the gap node structure at low energy. So, the very 
fact that the experimental observation cannot at present completely unravel the node structure of the gap function, 
reinforces the necessity of classifying all the possible symmetry breaking schemes. 

Moreover, the role of Fermi isotropic space beyond a zero-th order approximation in phenomcnological theories, is 
still supported by some recent studies: In HTS oxides, low symmetrical crystal fields (tetragonal and orthorombic 
as well as spontaneous strain, have weak infiuence on the temperature of the superconducting phase transition 
or on the penetration depth [Q. So, even if the spectrum of the interactions which are responsible for the pairing 



(whatever there origins may be) must be anisotropic in the crystal, it is worth analyzing the possibility that they act 
in a way not directly dependent on the crystal structure (like in isotropic space). 

The possible ground states of a high L superfluid has been the object of intense investi gati ons during the 60's and 
the 70's. ^From the solution of the state equations, Anderson and Morel ||l|] and Mermin |ll[| identified five different 
phases, through the minimization of a 4-th degree Landau potential. Schakel and Bias [ |12| analyzed the problem 
using only group theoretical arguments. Capel and Schakel ||l^, taking advantage of the results in ||l^, computed 
the possible ground states of condensates driven by p- waves. They also investigated the consequences of lowering 
the residual symmetry by means of strong spin-orbit interactions. Since, in this case, the Cooper pair is in a J = 2 
state [S ~1,L = 1), the order parameter is represented by a traceless symmetric tensor. The same holds true for an 
(5 = 0, L = 2) state order parameter, so their results may be directly applied to D-wave pairing. According to the 
authors of ref. [0|, eleven different phases are allowed. Their analysis, however, is questionable, since time reversal 
symmetry cannot be neglected in the theory of condensate states and the use of a 4-th degree polynomial Landau free 
energy strongly limits the number of phases that would be allowed by the symmetry of the system. 

Our aim is to give a definitive answer to the classification of possible symmetry breaking patterns in D-wave 
pairing Bose condensate, in the framework of the Landau theory of phase transitions. In this paper, wc shall mainly 
be concerned with the mathematical aspects of the problem, that will be solved making use of the geometrical 
invariant theory approach, proposed in ref. [Q. The strategy is to exploit a set of basic invariant polynomials of the 
symmetry group of the system, as fundamental variables in the description of the phase space of the system and in 
the minimization procedure of the free energy |15|| . In the realization of this program, we shall meet some substantial 
computational difficulties, that will be overcome by means of innovative procedures. 

II. THE GEOMETRICAL INVARIANT THEORY APPROACH TO SPONTANEOUS SYMMETRY 

BREAKING 

Let us briefly recall some basic elements of the geometrical invariant theory approach, or orbit space approach (see 
|16| and references therein) to the determination of possible patterns of spontaneous symmetry breaking. 

Let X € R" be a vector order parameter, transforming linearly and orthogonally under the compact real symmetry 
group G, and let <I>(a; x) be the G-invariant free energy, expressed in terms, also, of state variables a. The points 
xo{o:), where the function 4>a{x) = ^(a; x) takes on its absolute minimum, determine the stable phase of the system, 
whose residual symmetry is defined by the isotropy subgroup of G at xq, Gxq- Owing to its G-invariance, the free 
energy is a constant along each G-orbit, so, each of its stationary points is degenerate along the G-orbit through it. 
Since the isotropy subgroups of G at points of the same orbit are conjugate in G, only the conjugate class, [G^^q], of 
Gxg in G, i.e. the orbit-type (or symmetry) of the orbit through xq, is physically relevant. 

The set of all G-orbits, endowed with the quotient topology and differentiable structure, forms the orbit space, 
R^/G, of G and the subset of all the G-orbits with the same orbit-type forms a stratum of R"/G. Phase transitions 
take place when, by varying the values of the a's, the point a;o(Q;) is shifted to an orbit lying on a different stratum. 

If $(a; x) is a sufficiently general function of the a's, by varying these parameters, the point xoicx) can be shifted 
to any stratum of R"/G. So, the strata are in a one-to-one correspondence with the symmetry phases allowed by 
the G-invariance of the free energy. On the contrary, extra restrictions on the form of the free energy function, not 
coming from G-symmetry requirements (e.g., the assumption that the free energy is a polynomial of low degree), can 
limit the number of allowed phases for the system in its ground state. 

Being constant along each G-orbit, the free energy may be conveniently thought of as a function defined in the 
orbit space of G. This fact can be formalized using some basic results of invariant theory. In fact, the G-invariant 
polynomial functions separate the G-orbits, meaning that, for any two distinct G-orbits, there is at least a polynomial 
G-invariant function assuming different values on them. Moreover, every G-invariant polynomial can be built as a 
real polynomial function of a finite set, {pi{x), . . . ,pq(x)}, of basic polynomial invariants (integrity basis of the ring 
of G-invariant polynomials). Thus, the elements of an integrity basis can be conveniently used as coordinates of the 
orbit space points. They need not, for general compact groups, be algebraically independent, but can, and will, be 
chosen to be homogeneous polynomials in x. The number of algebraically independent elements in a minimal set of 
basic polynomial invariants is n — v, where v is the dimension of a generic (principal) orbit of G. Information on the 
number and degrees of a minimal set of basic invariants, and the degrees of the algebraic relations (syzygies) among 
them, can be inferred from the Molien function of G. 

Let us call qq the number of independent elements of the set {p}. The range of the orbit map, x i— > p{x) = 
{pi{x), . . . ,pq{x)) G R"^, yields a realization of the orbit space of the linear group G, as a connected semi-algebraic 
surface, i.e. a subset of R', determined by algebraic equations and inequalities. The orbit space of G will, therefore, be 
identified with a closed and connected region, S, of a go-dimensional algebraic surface, delimited by lower dimensional 



semi- algebraic surfaces. 

If an integrity basis has been determined, the equations and inequalities defining the orbit space of a compact group 
can be obtained from a simple recipe. It has been shown, in fact, that the orbit space of a reductive linear group can 
be identified with the semialgebraic variety formed by the points p € R', satisfying the following conditions i) and ii) 

i) p lies on the surface, Z, defined by the syzygies; 
ii) the q X q matrix P{p), defined by the relations 

n 

Pat{p{x)) = Y. ^^P-^""^ djPbix), \fx G R" (1) 

i=i 

is positive semi-definite and has rank go at p. 

The minimum of <I>(a; x) can be computed as a constrained minimum of the function, $(a;p), p g 5', defined by 

$(a;p(x)) = $(a;a;), Va; G R". (2) 

It has also been shown ||lj] that the state equation 

9$(a; x)/dxj =0, j = 1, . . . , n, 
determining the extremal points of <&(«; x), is equivalent to the following equation in orbit space: 

Y,Pabip)db^a;p) = 0, a^l,...,q, pes. (3) 

h=l 

Like all semi-algebraic sets, the orbit space of G presents a natural stratification. It can, in fact, be considered 
as the disjoint union of a finite number of semi-algebraic subsets of decreasing dimensions (geometrical strata), each 
geometrical stratum being in the border of a higher dimensional one, but for the highest dimensional stratum, 
which is unique {principal stratum). The geometrical strata are the connected components of the symmetry strata. 
The symmetries of two bordering strata are related by a group-subgroup relation and the orbit-type of the lower 
dimensional stratum is larger. 

In order to determine the minimum of $(a;p), p e 5, as a constrained minimum, one needs the relations defining 
the geometrical strata of R"/G'. These can be obtained from positivity and rank conditions on the matrix P{p) and 
from the syzygies, as already stated. 

The geometrical and symmetry strata are determined in the following way. Let W^'^'^ denote the (generally non 
connected) d dimensional algebraic subset of Z, defined by the relation 

W'^'^^ ^{pe Z\ rank(P(p)) = d, P{p) > 0}, (4) 

then, the geometrical d dimensional strata are the connected components of W^'^' and each d dimensional symmetry 
stratum, S'''''"-', a = 1,..., is the union of all the geometrical d dimensional strata, w'''^'"'^\ r — l,...ka, with 
the same orbit-type, [G*^'^''"-']. A representative in [G^'^'^'>] can be obtained as the isotropy subgroup at an arbitrarily 
chosen point x^'^'^'' of the image in R", VF^''''"-', of the geometric stratum W^'^'^\ The point x^'^'^'' can be obtained as 
a solution of the equation|j p{x) — p^''''''', where p^'^^^') is an arbitrarily chosen point in W^'^'^\ 

In the following, we shall classify and characterize all the allowed symmetry phases and possible phase transitions 
between contiguous phases, for I?-wave phase condensates in isotropic space. In particular, in Section III we shall 
identify the linear symmetry group, G, of these systems and a minimal set of basic polynomial invariants of G. In 
Section IV, we shall determine the geometrical features of the orbit space of G, i.e., its stratification (including 



^For p'"^ £ 5*, all the solutions, x, of the equation p{x) — p'"' form a unique G-orbit. The equation has no solutions for 

9(«) ^ S. 



connection properties and bordering relations of the strata) and the orbit-types of its strata. In section V, using 
an innovative method, we shall provide relatively simple expressions for the equalities and inequalities determining 
each stratum. These are essential for establishing the connection properties of the strata and make much easier the 
determination of the minima of a specific free energy function. A simple example of the determination of the minima 
of a general 4-th degree G-invariant polynomial and their stability against 6-th order perturbations, will be given in 
Section V. 

III. SYMMETRY OF THE ALLOWED D-WAVE CONDENSATE STATES IN ISOTROPIC SPACE 

The formation of D-wave condensate ground states breaks the symmetry of the isotropic S-dimensional space, which 
corresponds to the group O3 (g) Ui x (T), where O3 is the complete rotation group, Ui is the gauge group and (T) 
is the groupQ generated by the time reversal operator T. 

The symmetry of the allowed D-wave condensate ground states is defined by the relative values of the complex 
coefficients in the decomposition of the gap-function, A, in terms of spherical harmonics with L = 2: 



A(0,^)= }2 D^mYne,^). 

The set of functions {Y^,Y^*} yields a basis of a ten-dimensional (10 D) space, hosting a real representation of 
the symmetry group 03(g)Ui x (T). A general element, 7, of the group will be denoted by a triple 

7 = {P, e'^ e), 

where, p £ O3, < </> < 27r and e = —1, or +1 according as a time reflection is involved in the transformation, or not. 

In the following, we shall also use the shortened notations p for {p, 1, 1), Ui{(p) for (I3, Ui{(j)), 1) and T for (I3, 1, —1). 

The action of G can be transferred to a real irreducible action on the 10 D vector formed by the coefficients 

{D2, . . . , -D-2, D'^2j ■ • • 1 ^2}. The representation of G thus obtained can be realized in the lOD real vector space 

generated by a couple of two independent, real, second rank, symmetric, traceless tensors, X^, and X^, , i, j = 1, 2, 3, 
which can be considered as the real and imaginary parts of a complex 3x3 matrix ip, whose elements will be written 



in terms of five complex coordinates, z 



J- 



e R, (5) 



^-^ ^1 -^2 + ^ Z3 . (6) 
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The coordinates Da are connected to the Zj by the following relations: 

n izi + Z2 iz3 + Z4 iz3~ Z4 i Zi - Z2 , , 

The matrix ijj transforms in the following way under a general transformation, 7 — (p, (j), e) G G: 

7-^ = e'^pi^'pT, (8) 

where ip' — ^ or ^* , according as e = -1-1 or —1 and the apex T denotes transposition. As a consequence, the group 
G acts as a group of linear, real, orthogonal transformations on the vector order parameter x = (xi, . . . , xio) G R^". 
The kernel of the representation of G just defined is the group generated by the space refiection, so, it will not be 
restrictive to assume that the symmetry group is 



^Here and in the following we denote by {gi, . . . ,gn) the group generated by the set {gi, . . . ,gn}- 



G = SOa ® Ui X (T) 

and, when referring to G, in the following, we shall always mean this linear group acting in the vector space R^°. 

The group G has a trivial principal isotropy subgroup (the isotropy subgroup of generic points of R^°), thus the 
principal G-orbits have the same dimensions, four, as G and its orbit space, R^"/G, has dimensions qo — 10 — i — 6. 

The Molicn function of G, Mcirj), can be calculated in the form of an invariant Haar integral over G (see, for 
instance, |l^Jl^ ): 

where ^{g) is a normalized invariant measure on the group G, the integration is over the whole group G and g d G. 
In terms of the Euler angles 01, 9, 02 and of the gauge angle a, the integral can be written in the form of a sum over 
the two connected components of G 



-1 /'27r /'27r pTT 

167r^ Jo Jo Jo 



/ 1^ Jl [(i_^e*(fex+-))(^l_^e.(fex-a))] + J] [(l-r;e^'=x)(l + 77e^^-x)] 

■^" fc=-2 fc=-2 



(9) 



lVlG[ri) — 1^ .„2Vl _ Tl4^2/1 _ T^6^2/1 _ TiS'iS/i _ TilOVl — Til2A ' ^ ' 



where {e*'^'-^}_2<fc<2 are the distinct eigenvalues of the lOD rotation matrix and cos^ is a function of the Euler 
angles^: 

cos X = - [cos Q + COS(01 + 02) + COS Q COS(01 + 02) — 1] ■ (10) 

An explicit calculation of the integrals (see Appendix A) leads to 

a relation that can also be written in the following equivalent suggestive form: 

1 _ „16 _ „18 _ „20 _ 22 _ 24 _^ 26 _j_ „28 _|_ „30 ^ 32 ^ ^34 _ ^50 
(1 - r72)(l _ ^4)2(1 _ ^6)2(1 _ ^8)2(1 _ ^10)(1 _ ^12) 

Equations ( [ll| , [12| ) yield the following indications, whose validity has been checked through direct calculationsQ: 

1. A minimal integrity basis (IB) for the linear group G contains nine elements, {pi,...,pg} with degrees 
(di,...,d9) = (2, 4, 4, 6, 6, 8, 8, 10, 12). 

2. The invariants pi are connected by five independent syzygies of degrees 16, 18, 20, 22 and 24. The group is, 
therefore, non coregular. 

3. It is possible to find a non minimal 
homogeneous IB, {^i, . . . ,^6i ^i, • ■ ■ -Si], with degrees (2, 4, 4, 6, 6, 8, 8, 10, 12, 20), such that any invariant 
polynomial, F , can be written in the form [ pO| 

4 

^-Qo(0+EQ*(^)^- (13) 

i=l 

where the Q's, are polynomials in the algebraically independent polynomial invariants (^i, . . . ,^6)- 



^It has to be noted that, in principle, the denominators in (H), are polynomials in sinx and cosx; but, by they definition, 
they are also even functions of x, so, they can be expressed as polynomials only in cos x- 

*To justify the statements listed in the following items l)-3), starting from ( |ll| , |l2| ), one has to check that any G-invariant 
homogeneous polynomial of degree 20 in x G R^" can be written as a polynomial in pi (a;), . . . ,P'j{x). 



A possible choice for the IB's {p} and {^, 9} is the following: 

II I jS($3JJ2; PS = 5R Tr(^2^*)Tr(^3) (Tr, 



P2 



(14) 



and 

Ci^Pi, i = l,...,5, C6==P6+P7, /-^QN 

01^P6-P7, 02^ PS, 03^P9, di = {P6 - P7)P9- 

Thus, the most general non-equilibrium polynomial Landau potential, $(q!;p), can be written in the form 

^ = Qo + QlP7 + Q2PS + Q3P9 + Q'iP7P9- (16) 

Using the above definitions, the explicit form of the syzygies and of the P-matrix elements can be obtained through 
the following simple procedure. 

After defining the weight of an element, pi, of the IB to be equal to the degree, di, of the homogeneous polynomial 
Pi{x), and the degree of Yli p"' as ^^ Ui di, one writes down the most general homogeneous polynomial, Q^'^'{^, 9), 
of the relevant weight, d, in the weighted variables ^i, i — 1, ... ,6 and 9i, i — 1, . . . ,4. With the definitions 9o = 1, 
do = 0, the polynomial 6^''''(^, 9) can be written in the compact form 

3 6 

j=0 fc=l 

where, Nj = J2i=i "•* '^i ~ d-j, SnjA is a Kronecker symbol and the A's are real coefficient to be determined through 
the following conditions i) and ii): 

i) in the case of the syzygies S(i(Ci ^): d, — 16, 18, 20, 22, 24: The expression S£j(^(x), 6'(a;)) has to vanish identically 
in x; 

ii) in the case of the matrix elements Pa,b{£,, 9), a,b = 1, . . . , 9, d = da + di, — 2: The expression Pab{^{x), 9{x)) — 
X]i=i diPa{x) diPh{x) has to be an identity in x. 

The explicit expressions of the syzygies are reported in Appendix B and the elements of the matrix P{p) in Appendix 
C. In the following we shall always use the minimal IB {p}. 

On the algebraic surface, Z, defined by the syzygies in the 9 D space of the p's, the 9x9 symmetric matrix P{p) 
has rank 6, as it has to be, since the orbit space is six dimensional. 

IV. THE ORBIT SPACE OF THE LINEAR GROUP G 

As recalled in the Introduction, the orbit space of G can be identified with the range of the orbit map x i—> p{x), 
which is the same as the set of values of p G Z, that render the matrix P{p) positive semi-definite. 

The analytic conditions determining the geometrical strata, W^''"°''^\ are not difficult to specify, in principle, but, 
in the case of the higher dimensional strata, their explicit expressions (as sets of equalities and inequalities in the p's) 
become difficult to obtain, too large to be written down extensively and uneasy to handle, if they are derived only 
from rank and positivity conditions of the matrix P(j>). An elegant way to bypass this difficulty is the following. 

One starts with the relatively easy determination of the geometrical strata of dimensions < 2 {i.e., the connected 
components of the sets W'^'^\ d = 1,2, defined in (Q)), starting from the positivity and rank properties of the matrix 
P{p) and from the syzygies. The subsequent step consists in the identification of the orbit-type of each geometrical 
stratum. This allows to identify the connected components, W'''^'°''^\ r = 1, . . ., with the same orbit-type, [G'-'*'"''], 



and belonging, therefore, to the same symmetry stratum^, S^'^'"\ All the groups G^^'"^ turn out to be finite groups 
of low order (< 8). 

The determination of the higher dimensional strata bordered by a given 2D stratum, 5"*^^'"^ can be obtained 
starting from the identification of their orbit-types, which are necessarily contained in the orbit-type of S^'^'"\ The 
task is easily accomplished by selecting, out of the set of maximal subgroups of each G^^'"^ those which are isotropy 
subgroups of G. To this end, the following criterion can be used. Let H be one of the maximal subgroups of a G^'^'"-' 
and denote by Vh the linear subspace formed by all the points x E R^", which are stable under H. Then, H is an 
isotropy subgroup of G, if and only if the following condition is satisfied: 

Maxrank(P(a;)) < d. 

If this condition is satisfied, when x spans Vff , the point p(.t) e R^ spans the topological closure of a stratum 5^'' '" •* 
and 

dim(S'(''''"')) = Max rank(P(a;)). 

To get the full set of higher dimensional strata, the last part of the procedure just described has to be repeated for 
each newly determined stratum. 

The derivation of the equalities and inequalities determining each of the higher dimensional strata requires a more 
sophisticated method pi| , which we shall explain and exploit in the following section. 

Before resuming the results obtained with the realization of the program just presented, let us note that, owing to 
the linearity of G, its isotropy subgroups coincide at the points x G R"'^^ and kx, 7^ k G R. Thus, taking also into 
account the homogeneity of the polynomials Pi(a;), the points (pi, . . . ,_P9) and (I,p2/Pi^ , • • ■ iPg/Pi^ ) lie on the 
same stratum and the equalities and inequalities defining the strata can be written as homogeneous relations in the 
weighted variables pi. It will be useful, therefore, to introduce the definitions 

Vi^VilV^X , P= (P2,...,P9). (17) 

The relations determining the ID and 2D strata can be easily obtained by selecting those solutions, p, of the 
condition rank(P(p)) — 1 and, respectively, rank(P(p)) = 2, at which the matrix P(p) is positive semi-definite. In 
this way, one obtains five connected components, VF(1'^), A — 1, ... ,5, of W'^^'^ and seven connected components, 
^(2,-4)^ A = 1+, 1_, 2, 3+,3_,4, 5, of VF*^^-'. The parametric equations defining this semi-algebraic sets are listed in 
Tables | and g. _ 

In order to identify the distinct ID and 2D symmetry strata, for each connected component W^'^'^> ^ d = 1,2, we 
have picked up a point ■p^'^^^^ and, for each p('*''^), a "simple" solution, x^*^'"^-*, of the equation p[x) — p^'^'^\ Then, 
at each point a;^''''^) we have determined, and compared for conjugation, the isotropy subgroup of G, G^"^'"^' . These 
subgroups turn out to be non conjugate in G, but for G*-^'^+-', which is conjugate to G'^^'^^\ and G'^^''^+\ which is 
conjugate to G'^^'^^-'. Thus, there are five distinct ID and five 2D symmetry strata, S^'^'°'\ with orbit-types [G'^'''"-'], 
d = 1, 2, a = 1, . . . , 5. 

After choosing an element, H, in each class [G*^'^'"-'], it is easy to determine the subspace Vh C R^°, formed by the 
vectors which are if -invariant. 



The results we have obtained are listed in Tables III and IV and the bordering relations among the strata are 
illustrated in Figure H. The definitions of the rotation matrices appearing in the tables and in the following text are 
recalled in Appendix D. We have used standard notations for the corresponding geometrical transformations (see, for 
instance, p^). 



^There is a unique zero dimensional stratum, corresponding to the point p = p(0) = 0, whose symmetry is [G]. It will be 
ignored in the following. 



V. RELATIONS DEFINING THE STRATA WITH DIMENSIONS > 3 

In order to derive the relations defining the higher dimensional strata we shall exploit the following results [ pi[ , 
that we shall briefly recall, without proof. Let H be an isotropy subgroup of G, S[h\ E^nd S[h\ the associated strata 
in R^" and, respectively, Ti}^ /G and let Stab(i/, G) be the stabilizer of H in G: 

Stah{H,G) ^ {g e G \ g H g-^ ^ H}. 
We have denoted by Vh the linear subspace formed by the iJ- invariant vectors of R^°: 

Vh ^{xeR^" \hx = x, yhe H}, (18) 

and we shall denote by V^ the subset of Vh formed by the points at which the isotropy subgroup of G is H: 

V^"^ ^{xe Ri" I G, = H}. (19) 

The set V^ ' is the intersection of Vh with the stratum of orbit- type [H] . 

The group Stab(iJ, G) acts linearly on Vh', let us call H the linear group defined by this action. The i/-orbit through 
X S Vh is the intersection of the G-orbit through x with Vh- There is, therefore, a one-to-one correspondence between 
the G-orbits lying in the closure of S[h] Sind the H orbits. This correspondence can be restricted to a one-to-one 
correspondence between the interior points of the closure of S^h] and the points of the principal stratum of Vh/H. If 
S[H] is connected, the one-to-one correspondence is between the principal stratum of Vh/H and S[h], otherwise the 
inverse image of S[h] in the correspondence, reduces to the points of the principal stratum of Vh/H, corresponding 
to if -orbits through points of V^ '. Thus an IB, {A}, of H can be used to parameterize the points of Syn]'- 

The range of the parameters A will be determined by the positivity conditions of the P-matrix associated to the IB 
{A}, P'^'(A), and convenient additional conditions obtained from positivity and rank conditions of P{p), if S[h] is 

not connected. To this end it will be worth noting that, for p e Sih] ■ the matrix P{p) can be written in the following 
form: 

P{P{X)) = J^{X)P("\X)J{X), (20) 

where, J(A) is the Jacobian matrix 

..W^^. (21, 



Owing to a well known theorem in matrix theory, from ( |20| ) one obtains the following upper limit: 

Tiink(p{p{X))] < rank(J(A)). (22) 

Below, we shall use this approach to derive rational parametric equations for all the strata of S with dimensions 
three and four. The results obtained in this way, starting from each of the strata of dimensions > 2, are resumed in 



Table IV and in the lower part of Table III 



1. Stratum 5'^'^' 

The stabilizer in G of the group H = G^^.i) = (Ga^, C2^^T), is the group Stab(i7) = (Ga^, C2a, Ui{t:),T). 

An easy calculation shows that the space Vh is defined by the equations Xi = 0, i = 1, 3, 4, 7, 8, 9, 10. The 
elements of H act on the Xi, i = 2, 5, 6 by changing the signs of one or two or all of these coordinates, so, an IB of H 
can be chosen to be 



Ai=a;5, A2 = a;2, A3 = Xg (23) 

and the expression of the p's in terms of the A's turns out to be the foUowing: 

Pi — Ai + A2 + A3, 

P2 = (3AH6AiA2+3A2 + 2AiA3 + 18A2A3+3Ai)/6, 

P3 = (Ai + A2 - x^y, 

P4 = Ai(Ai-3A2 + 3A3)V6, 

P5 = Ai(Ai-3A2-A3)V6, (24) 

P6 = Ai (Ai - 3A2 - A3) (Ai + A2 

P7 = Ai(Ai-3A2-A3)2(Ai+A2 

Ps = Ai(Ai ~3A2- A3)(Ai + A2 

P9 = Ai(Ai+A2-A3)3(Ai-3A2 + 3A3)V6, 

with 

Ai, A2, A3 > 0. (25) 



A3) (Ai - 


-3A2 


+ 3A3)/6, 


-A3)/6, 






A3)' (Ai 


-3A 


2+3A3)/6, 



For values of the A's satisfying (Eq), the values of the p's defined in (24) render P{p) everywhere positive semi-definite 
and of rank three. We can conclude, therefore, that they yield parametric equations for the 3D stratum S^^'^'. Being 
the continuous image of a connected set, the stratum is connected. 

By eliminating the parameters A's from (p4 25), one obtains the following relations: 



P5 = (2pi^ -2p2~P'i) {2 pi'^ - 6 p2 ~ 2 pi ^/p^ ~ ps) T, 

P6 = -'Sy^ {2pi^ ~2p2~P3) {2pi^ -6p2~2piy/p^~p3) {2pi^ ~2p2 + 2piy/p^ + p3) T, 

P7 = -V^ (2pi' - 2P2-P3) (2pi' -6p2-2pi^-p3) r, 

Ps = 3p3 {2pi'^ -2P2-P3) {2pi^ -6p2-2pi^-p3) {2pi^ -2p2 + 2pi^ + P3) T, 

P9 = -9Vp! (2pi2-2p2-P3) (2pi2-2p2 + 2piv^ + p3)V, 



< P3 < pI 

2pl+p3 <6p2 < 3{2pi-p3), 



,2 I ^„ / R^. / q /'o^2 ^\ (27) 



where 



r-i - 64 (pi + VP;)^ (28) 



and the square root of P3 has to be intended in an algebraic sense. 

2. Stratum S<^'^' 

The stabilizer in G of the group H — G^^-^' = {C'2x,C2z), is a group of 00 order: Stab(iJ) — {Cix, Caz,T) x Ui. 

The space Vh is defined by the equations xi = 0, i = 1, 3, 4, 6, 8, 9. The transformation properties of the 
coordinates x^, i = 2, 5, 7, 10, under the transformations of H can be conveniently described in the following way. 
Let us define: 

r]i^ X2 + ix^ + i{xT + ixiQ), r]2 = X2 - i X5 + i{x7 - ixio), (29) 

then 

i) Under a gauge transformation, Ui(0): (?7i,?72) -^ (e"^»7i, e*'^?72); 
ii) under time reversal T: (771,772) —* ('72i'7i); 
iii) under a transformation C42:: (771,772) -^ (e''^^'^ 772, e"*'^/'^ 771); 
iv) under a transformation C4z: (771,772) — > (—772,-771). 



The group H is coregular and admits the foUowing minimal IB: 

A2 = {\m? -\m?Y/'^^ = {xhX7~x2Xiof, ,„„^ 

A3 = 5R[(r7i%*)'] - {x2^ - x^ + x,^ ~ x^^) {x2^~Ux^x^+x,^ + 2x2^x,^~ ^^^^ 

2 ^5^ xy"^ + xy** - 24 a;2 ^5 a;? xw - 2 a;2^ xio^ + 2 xs^ xiq^ - 14 xy^ xio^ + xw^) . 

It has to be noted that the invariant (5 [(771 772)'^] ) can be expressed as a polynomial in Ai, A2, A3. 
The corresponding P-matrix turns out to be 

4Ai 8A2 I2A3 

8A2 4A1A2 I (31) 

12 A3 36Ai(A?-4A2)^ 

and the p's can be expressed in terms of the A's in the following form: 

Pi = Ai, 

P2 = (3Ai2-4A2)/6, 



P3 = Ai^-4A 



2, 



Pi = (Ai3 + 12AiA2-A3)/12, 

P5 = (Ai3-4AiA2-A3)/12, (32) 

P6 - (Ai^-16A2'-AiA3)/12, 

P7 - (Ai^-8Ai2A2 + 16A2'-AiA3)/12, 

PS = (Ai^-8Ai3A2 + 16AiA2'-Ai2A3-4A2A3)/12, 

P9 = (Ai*^-12Ai''A2 + 48Ai^A2^-64A2^-Ai^A3-12AiA2A3)/12, 

where the parameters A have to satisfy the conditions: 

Ai >0, 

< 4A2 < A?, (33) 

Ai < (A2-4A2)'. 



For values of the A's satisfying (P3|), the values of the p's defined in (32) render P{p) everywhere positive semi-definite 
and of rank three. We can conclude, therefore, that they yield parametric equations for the 3D stratum S'^^'-'^^ Being 
the continuous image of a connected set, the stratum is connected. 

By eliminating the parameters A from (|32| , p3| ) one obtains the following relations: 

P2 = f2pi2+p3)/6, 

P5 = (-Pi^ +PlP3 + 3p4) /3, 

P6 = {-Api'^ + hpi'^pz-p3^ + I2pipi) /12, , . 

P7 = -4pi4 + 3pi2p3+p32 + 12pip4)/12, ^^^^ 

PS = (-4pi'' + 5pi^P3-piP32 + 12pi2p4-6p3P4) /6, 

P9 = -16pi6 + 24pi4p3-9pi2p32+p33 + 48pi3p^_36p^P3P ^/12 



and 



pI< [l2p4-(4p?-3p3)Pi] . ^ ^ 



3. Stratum ^f^'^' 

The stabilizer in G of the group H = G^^'^) = {C2x,T), is the group Stab(iJ, G) = Of x (C/i(7r),T). 
The space Vf/ is defined by the equations Xi — O.i = 1, 4, 7, 8, 10. The transformation properties of the coordinates, 
Xi, i = 2, 3, 5, 6, 9 of Vu are easily derived from (|^) and can be put in the following advantageous form. Let us define 

V3X2+X5 ..X2-V3X5 . 

m^ 2 ' V2^X3+l , T]3=Xe-lXg, (36) 

then, 

10 



i) under a proper rotation, 02(0): (771,772, 773) ^ (?yi, e ^'"^ 772, 6**^773); 

ii) under a gauge transformation, Ui(7r): (771, 7^2, 773) -^ (-7/1,-772,-773); 
iii) under a transformation C2z'- ('yi, '72,%) — > ('7i, — '^2: 'ys); 
iv) under time reversal T: (771, 772, 7/3) — > (?7i, 772, -773)- 
The group iJ is not coregular and admits the foUowing minimal IB: 

Ai = 3|77i|2 = 1(73x2 +2:5)', 

" +4x1-2 V3x2a;5 +3x| 



A2 — |??2| — 4 1^2 T" '±■^3 ~ ^ v«J J-2 -i-s -r 0x5; , 

A3 = |??3p ^xe^ +X9'^ 

A4 = %/3 77i3(772 77|) = ^ (73x2 +a;5) (a;2 xg^ - n/Sxs xg^ - 4x3 xe xg- 



(37) 



X2X9' 



V3 



X5X9 



[3 (7/277!)] = J (X2X6^ - a/3x5X6^ -4X3X6X9 - X2X9^ + VS X5 Xg^ 



with the syzygy 



and the P-matrix 



/12Ai 




6A4 
V 





4A2 









4 A3 



A4 — A1A5 



6A4 
2A4 
4A4 



0. 



(38) 





4 As 

8 As 



2A4 4A4 4Ai A2A3 + Ai A3^ + 2A3(4A2 + A3)A4 

4A5 8 As 2A3(4A2 + A3) A4 4 A3 (4 A2 + A3) As / 



The conditions P(A) > and rank[p(A)') 
for the A's: 



(39) 

4, added to (BSh, yield the following restrictions on the acceptable range 
Ai, A2, A3 > 0, 



< As < A2A2, 
Ai + As > 0. 

The p's can be expressed in terms of the A's in the following form: 

Pi - (Ai+3A2 + 3A3)/3, 

P2 = (Ai^ + GAi A2 + 9A2^ + 14AiA3 + 18A2A3 + 9A3^ + 24A4)/18 

P3 - (Ai+3A2-3A3)V9, 

P4 = (4Ai^-72Ai^A2 + 324AiA2^-36Ai2A3 + 324AiA2A3+ 



(40) 



P5 
P6 
P7 

Ps 

P9 



\i — 72 Ai A2 
2 



324 Ai A2^ + 12 Ai^ A3 - 108 Ai A2 A3+ 



81 Ai A3^ + 108 Ai A4 - 972 A2 A4 - 486 A3 A4 + 729 As) /648, 

(4A3 

9Ai 

(Ai 

IO8A1A2A3-27A1A3 



9 Ai A3^ - 36 Ai A4 + 324 A2 A4 - 54 A3 A4 + 81 Ag) /648, 
(A1+3A2 



3 A3) (4 Ai^ - 72 Ai^ A2 + 324 Ai A2^ - 12 Ai^ A3 + 

^ + 36 Ai A4 - 324 A2 A4 + 162 A3 A4 - 243 As) /1944, 
(Ai + 3 A2 - 3 A3) (4 Ai^ - 72 Ai^ A2 + 324 Ai X2^ + 12 Ai^ A3- 
108 Ai A2 A3 + 9 Ai A3^ - 36 Ai A4 + 324 A2 A4 - 54 A3 A4 + 81 As) /1944, 
(Ai + 3 A2 - 3 Aa)^ (4 Ai^ - 72 Ai^ A2 + 324 Ai A2^ - 12 Ai^ A3+ 
108 Ai A2 A3 - 27 Ai A3^ + 36 Ai A4 - 324 A2 A4 + 162 A3 A4 - 243 As) /5832, 



(41) 



(A1+3A2 



3 A3)' (4Ai' 



72Ai^A2 + 324AiA2^ 



36Ai A3 



324 Ai A2 A3 + 81 Ai A3^ + 108 Ai A4 - 972 A2 A4 - 486 A3 A4 + 729 As) /17496. 



For values of the A's satisfying (|38| ) and (40), the values of the p's defined in ( [4l[ ) render P{p) everywhere positive 
semi-definite and of rank four. We can conclude, therefore, that they yield parametric equations for the 4D stratum 
S''-'*'^'. Being the continuous image of a connected set, the stratum is connected. 
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For Ai 7^ [A5 ^ 0], the syzygy can be solved with respect to A5 [Ai] and the expression one obtains can be 
substituted into (^, so, in this region, (Ai, . . . , A4) [(A2, . . . , A5)] play the role of local coordinates for the manifold 
underlying the stratum. Since, owing to (pO|), Ai and A5 cannot vanish simultaneouly, the whole stratum is covered 
by two rational charts of local coordinates. 

The elimination of the parameters A from (BllBn) leads to cumbersome relations, which it is not worth writing 
down. 

4. Stratum S*'*'^' 

The stabilizer in G of the isotropy subgroup H = G*^'*'^' = ((^22), is the group Stab(_ff) = O2 x Ui. 
The space Vr is defined by the equations Xi = 0, i = 3, 4, 8, 9. The transformation properties of the coordinates, 
Xi, 1 = 1,2, 5, 6, 7, 10 of Vh can be put in the following advantageous form. Let us define 

x-j + ixii + i{x2 + ixi) x-j + ixQ- i{x2+ixi) , ^^. 

'ni = x^ + zxio, m = -= , 773 = -= , (42) 

then, 

i) under a proper rotation, 02(0): (?7i,?72, Jys) ^ (?yi, e^^*'^, e^^*"^ 773); 
ii) under a gauge transformation, Ui((/i): (?7i,?72,%) -^ (e**^ ?7i, e~*"* 772, e*'^7y3); 
iii) under a transformation C2x'- ('7i,??2,??3) —^ ('/i : ^73 1 ^2 ) 5 
iv) under time reversal T: (771, 772, %) — > (?7i, -%, -?72)- 
A minimal IB of iJ, which is coregular, can be chosen to be 

Ai = Y.]=i \V]? ^ xi^ + ^2^ + x^"^ + X(i^ + ^7'^ + xio^ ^ 
A2 = \vi? ^xl+xjo, 

A3 = {\rn\'~\V3\')y^={x2Xe-x^xr)\ (43) 

A4 = 2 3? (772 772773) ^ -xlxl-xlxl+xlxl+xlxj- 
4: Xi X5 Xq Xio -4X2X5 X-j XiQ + x\ x\q + xl x\q - xl xIq - Xj xIq. 

The P matrix relative to the IB {A} has the following form: 

/ Ai A2 2 A3 2A4 \ 

V2A4 A4 A2(A2-AiA2-4A3)/ 

The conditions P(A) > and rank! P(p) j = 4 yield the following restrictions on the acceptable values of the A's: 

< A2 < Ai, 

< 4A3 < (A1-A2)', (45) 



A4 < A2 



\2 



(A1-A2) -4A3 . 

The expressions of the p's in terms of the A's turn out to be the following: 
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Pi 

P2 
P3 
Pi 
Pb 



P7 



Al, 

(3Ai 
Ai^- 
A2 



P9 



- 2 Al A2 + 2 A2 
2 Al A2 + 2 A2 - 



- + I2A3-2A4 
4A3-2A4, 



/6, 



(9Ai 



18AiA2 + 10A2^-36A 



= (5Ai A2-I4A1A 

>2 



= (6 



Al A2 



-ITAi^Aa" 



o Al A4 -f- 2 Al A2 A4 
(4Ai3A2 



llAi^Aa" 



5Ai A4 + 10AiA2A4-4A2^A 



(3Ai4A2 + 2Ai'A2' 



> Al A2 A3 



+ 80 Ao^ A3 + 



8 A2^ A4 + 24 Al A3 A4 - 16 A2 A3 A4 - 4 Ai A4^ 
(33Ai'^A2^-132Ai= 



6A4)/6, 

10 A2^ - 4 A2 A3 - 4 Al A4 + 6 A2 A4) /6, 
14AiA2^-2A2''-24AiA2A3 
12A3A4-2A42)/6, 
-2A2''-16AiA2A3 + 28A2^A; 
4A3A4+2A42)/6, 
+ 52 Al A2'^ - 20 A2^ 
IBAi^A 
-12A2A4^)/6, 

5 



4A2'A4 
6 Al A2 



36Ai2A2^ 



) A2 A3 - — b Al A4 
U-4. 

2 X 4 



20A2"A3 



24Ai"A2A3- 
2 A4 — Id Al A2 A4- 



A2'' + 168Ai^A2^-72AiA2" 



528 Al A2-' A3 - 144 A2^ A3 



528A2^A3^ 



4A2®-264Ai^A2^A3+ 
9 Ai"^ A4 + 36 A? A2 A4 - 60 \l \l A4- 



i Al A2^ A4 - 12 A2* A4 + 72 Al A3 A4 - 144 Ai A2 A3 A4 + 96 A2-' A3 A, 



144 As-" A4 



12Ai^A4^ 



24Ai A2 A4 



4— 
12 A2^ A4^ + 48 A3 A4^ - 4 A4^) /6. 



(46) 



In the range of values for the A's defined in (|45|), the values of the p's defined in ( |46[ ) render P{p) everywhere 
positive semi-definite and of rank four. We can conclude, therefore, that they yield parametric equations for the 4 D 
stratum S^*'^'. Being the continuous image of a connected set, the stratum is connected. 

The elimination of the parameters A from (Hq,B3) leads to cumbersome relations, which it is not worthwhile writing 
down. 



5. Principal stratum S^'''^' 



Implicit equations for the principal stratum are yielded in R^ by the syzygies relating the invariants pi, i = 1, . . . , 9, 



which can be easily obtained from the syzygies (B1-B5), using ( [l5| . At the points of 5p — S'^^'^^ the matrix P{p) has 
to be positive semi-definite and its rank has to be 6. 

The method used for the singular strata does not help to obtain parametric equations for Sp and to attain the goal 
one has to resort to less elegant and more pedantic procedures. A possibility is to use the invariants of an IB of the 
subgroup SO3 of G as parameters in terms of which to express the p's, but in this case, great care has to be paid to 
the determination of the range of the parameters that render one-to-one the correspondence with the points p (£ Sp. 

The Molien function for the subgroup SO3 of G can be obtained as a byproduct of the calculation of the Molien 
function for G (see Appendix A): 



MsoM 



l+rj^ + r]^ 



which, after multiplying numerator and denominator by (l — 77'*), becomes 



MsoAv) 



l-ryi 



(1-7^2)^ (l-ryS)^ (1-7?^) 



(47) 



(48) 



Equations (B^,E8[) indicate that a minimal IB for the subgroup SO3 of G is formed by eight invariants, with degrees 
(2, 2, 2, 3, 3, 3, 3, 4) related by only one syzygy of degre 12. A possible choice is the following: 



Al 
A2 
A3 
A4 



Tr [tjjip*) , 



^ 


"Tr (^-2^1 


3 


Tr W\ 


5R 


Tr \i>^) 



5 


^ 5 


Tr W)] , 


6 


= 5R 


Tr (7/>2 7/>* 


7 


== 5 


Tr (7/;2 i,* 


8 


= Tr 


(ipip* ipip 



(49) 



Only the first seven elements of the basi s ar e algebraically independent. The explicit expression of the syzygy relating 
the A's is reported in Appendix B (see (B7)). 

The P-matrix relative to the IB {Ai, . . . , As} is is the following and for values of the A's satisfying the syzygy 
it has rank seven: 
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Pi, = 2 df"^ K, d[^^ = (2, 2, 2, 3, 3, 3, 3, 4), i = 1, . . . , 8, 

P22 = 4Ai, P45 = 0, 

P23 = 0, -P46 = i(2^i^2 + A2'-A3'), 

P24 = 6A6, Pi7 = -(Ai-A2)A3, 

P25 = 6A7, hs = 2(2AiA4 + A2A6-A3A7), 

P26 = 2(A4 + 2A6), P55 = I (6Ai2 - A2' - A32 - GAg) , 

P27 == 2(A5-2A7), ^56 = (Ai + A2)A3, 

P28 = 4A1A2, i^57 = H2A1A2-A22 + A32), (50) 

P33 = 4Ai, P58 - 2(2AiA5 + A3A6 + A2A7), 

P34 = -6A7, Pee = j^ (2Ai' + 8A1A2 + 5A2' + 5A3' + eAg) , 

P35 = SAe, Pe? = 3A1A3, 

P36 = 2(A5 + 2A7), Pes = -| (3A2A4 + 3A3A5 - UAiAg + 2A2A6 + 2A3A7) , 

P37 = -2(A4-2A6), P77 = j^ (2Ai2 - 8A1A2 + 5A22 + 5A3' + GAg) , 

P38 = 4A1A3, -P78 = f (3A3A4 - 3A2A5 - 2A3A6 + I4A1A7 + 2A2A7) , 

P44 = f (6A? - A2 - Ai - 6Ag) , i^88 = -| (SAi^ - 2A4' - 2A5' + 2A6' + 2A7' - 9AiAg) . 

The connection between the p's and the A's can be immediately obtained from their very definitions: 

Pi = Ai, 
P2 = Ag, 
P3 = A2+A2, 

Pi = xj + xl 

P5 = Ai + A?, (51) 

Pq = A2 A4 Ag + A3 A5 Ag — A3 A4 A7 + A2 A5 A7, 

P7 = A2 Ag +2 A3 Ag A7 — A2 A7 , 

Pa = A2 A4 Ag — A3 A4 Ag + 2 A2 A3 A5 Ag + 2 A2 A3 A4 A7 — A2 A5 A7 + A3 A5 A7, 

P9 = A2 A4 — 3 A2 A3 A4 +6 A2 A3 A4 A5 — 2 A3 A4 A5 — A2 A5 +3 A2 A3 A5 . 

As aheady noted, there is not one-to-one correspondence between the points A = (Ai, . . . , Ag), lying on the surface 



determined by the syzygy (B7) and rendering positive semi-definite and of rank seven the P- matrix associated to the 
integrity basis {A}, and the points p of the principal stratum of R^°/G. In fact, the action of the group G in R^" 
induces a linear action of G on the variables A^ . With the definitions 

77i=A2+iA3, ??2 = A4 + iA5, % = Ag -I- i A7, (52) 

the variables A and 77 have the following transformation properties: 

i) Ai and Ag are invariant under the whole group; 

ii) under gauge transformations, Ui(0): (?7i, ?72, %) -^ (e^**^ ?7i, e^"^?72, e*"* 773); 

iii) under time reversal T: (?7i,?/2, ^73) —^ (^1 5^21^3 )• 

The one-to-one correspondence problem can be solved by finding a criterion enabling to select a unique point in 
each orbit of the action of G in the A space. To this end, let us meke to preliminary remarks; 

i) A direct check shows that for rji ~ and general values of 772 and 773 the rank of P{p{X)) is six. 

ii) Using ( p2|) one easily realizes that rank! P(p(A)) ) <rank(J(A)) < 6, for A3 = A5 = A7 = 0. 

Now, let us choose any orbit, J), in the principal stratum of R^^/G, i.e., ap & Sp, and let A be such that p{X) = p. 
At least one of the parameters fji, i = 1,2,3, has to be different from zero. Let us say that fji, is the first one, in 
lexicographic orderH. Then, on 0, there will be a point A', such that 77^ is real and positive; this point can be attained 



®In a lexicographic ordering, {vi, . . . ,Vk) > {v[ , . . . , v'f,) means that the first non vanishing component Vi — v'^, i = 1, . . . , fc, is 
positive. 
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from A by means of a convenient gauge transformation. Since at least a fjj, j — 1,2,3, has to be ^ 0, this means 
i < 3. At this point we still have at our disposal time reversal transformations to fix the sign of the imaginary part 
of one of the fji, i > i (recall that at least one of the 3(%) has to be ^ 0). So, if j is such that '^{fjj) is the first (in 
lexicographic order) non vanishing imaginary part of the 77's, we can choose 3(^j) > 0. 

After introducing a lexicographic order in the real vector spaces generated by the vectors formed, respectively, with 
the real and imaginary parts of the 77's, the criterion that wc have devised can be resumed in the following simple 
additional conditions on these vectors: 

A3 = 0, (A2, A4, Ag) > 0, (A3, As, A7) > 0, and A5 = for A2 = 0. (53) 

Using these conditions, the expressions { pV{) of the p's can be simplified as follows: 

Pe = A2 (A4 Xq + A5 A7) , 

Pj = A2 (A| - Af) , . 

PS = Ai(A4A6-A5A7), ^^^^ 

P9 = A3(Al-Ai). 



Pi 


- Ai, 


P2 


= Ag, 


P3 


= Ai, 


Pi 


= Al + Ai, 


P5 


= Ai + A?, 



where the A's are required to satisfy the syzygy (B7), the conditions rendering the matrix P{p) positive semi-definite 
and of rank six and (A5, A7) > 0. 

VI. AN EXAMPLE: MINIMA OF A GENERAL FOURTH DEGREE LANDAU POLYNOMIAL 

As a simple example, we have also calculated the minimum of a general fourth degree polynomial free energy 

3 

' 2 , ' 



^*'Hp)-^P? + E"jPJ' a^eR; i = 0,...,3, (55) 

in the additional assumptions that it is bounded below and has a local maximum at the origin (ai < 0). In ( p5| ) the 
a's are phenomenological parameters. 

Recalling the definition pi = pi/pi , the polynomial $('*'(p) can be put in the following form: 

2 
$(4)(p)^^A + a,p,, (56) 



2 



where we have defined 



A = do + 2a2P2 + 2a3P3. (57) 

Since, owing to its definition (|l4|), pi ranges over the whole non negative real numbers, in the orbit space, the 
polynomial <I>(^^ (p) is bounded below (in the assumption ai < 0) if and only if the minimum of A, (5, evaluated in the 
section pi = 1 of the orbit space, is positive. Being the minimum of the r.h.s. of (|56|), thought of as a function only 
of pi > 0, equal to — ai/(2A), the absolute minimum of $(*^(p)| „ is —a\/{26). In this way, the determination of 

the minimum of $'^'(p) is reduced to the calculation of 5. 

The absolute minimum of 5 in each singular stratum can be easily computed using the equations of the strata. For 
the principal stratum the determination of the minimum of $'■^•'(73) as a constrained minimum is difficult and it 

is easier to solve equation (g) and to check subsequently if the solutions lie in the principal stratum. A comparison 
of the values of the minima in the different strata, in order to determine the absolute minimum, leads to the results 
resumed in Table |V| and illustrated in Figure g, where the denomination of the strata introduced in the text has been 
used. Owing to the low degree of the polynomial defining $'^'(p) in (pa), and the consequent low number of free 
parameters a, the absolute minimum exhibits strong degeneracy, particularly for special values of the a's. If these 
special values are excluded, spontaneous breaking of the symmetry can generate only five distinct phases out of fifteen 
permitted in the G-symmetry; some of them are unstable. For no non trivial values of (a2,a3) does the absolute 
minimum lie on the stratum S^'^'^' . 

For general values of the a's, our results are in agreement with Mermin's ones [0. Let us add a few words about 
the perturbative stability of the three degenerate phases in the region TZ2 (see Table ^. For (ao, a2, a3) S TZ2, the 
addition to the free energy, $(''\ of a "small" perturbation, consisting in an invariant polynomial of degree six: 
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e(6)=a4P4 + a5P5, (58) 

splits the three degenerate minima determined by the 4-th degreee termQ This is easy to check, at least in the 
additional assumption that the perturbation leaves the absolute minimum in one of the strata corresponding to the 
degenerate phases. In fact, at the first perturbative order, one obtains from Tables | and || the following shifts, Q/J,.), 

in the values of the 6-th order free energy at the points where $'■'*•' (p) takes on its degenerate absolute minimum under 
consideration: 

(6) a? aj + as „(6) „ p.(6) a? . , ^ , ,.o\ 

Q(i.2) = ^ 6 ' ®(i.3)='^' ^(2,4) = ^ ("4 + a5)C, (59) 

where < ^ < 1/6. 

Since —ai/6 > 0, the absolute minimmn will be perturbatively stable on S**-^'^^ or, respectively, 5*^^'"^^ according as 
a4 -I- as is negative or positive. 

As stressed in the Introduction, the difficulties mentioned above can be overcome if one puts less restrictive upper 
limits to the degree of the polynomial describing the free energy. It is trivial, for instance, to realize that the following 
class of bounded below polynomial functions pq ] have a vanishing maximum at the origin of R^'' and display an 
absolute minimum at the arbitrarily chosen point p (z S: 
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[(P* -Pi)^"' --P?"' : (60) 



where the a's are positive constants and the n's are positive integers. 

The physical implications of our results and the derivation of a more realistic form of the free energy will be 
discussed in forthcoming papers. We shall limit ourselves, here, to note that, the D-wave content in terms, e.g., of d- 



or extended s- waves (mentioned in the introduction) can be directly read off the fourth column of Table III , for each 
of the possible symmetry phases of the system under consideration. 

VII. CONCLUSIONS 

Through a detailed and rigorous analysis, we have determined and characterized the possible ground states of D- 
wave condensates in isotropic space, which are rilevant both in high Tc superconductivity and in ^He phase transitions. 
The problem has been formulated and solved in the framework of geometric invariant theory, that, as noted long ago 
p4| , is the natural mathematical setting in which to deal with minimization of invariant potentials, which are plagued 
by degeneracy problems. This enabled us to complete and/or correct results previously obtained by other authors, 
using different methods. 

Until now, the kind of approach we have used had only been applied to relatively simple models, with coregular 
symmetry groups, in condensed matter p3| and in elementary particle physics [E4|. In this paper, we have fully 
characterized the geometry of the orbit space of the linear group SOa^Ui x (T), acting in the space R^", a highly 
non trivial group, with a non coregular minimal integrity basis, formed by nine elements with degrees up to 12. For 
this group we have determined the explicit form of the elements of an integrity basis and of the syzygies relating 
these elements, the equations and inequalities determining the orbit space and its stratification and the correspond- 
ing isotropy subgroups chains. The minima of a general fourth degree polynomial invariant free energy have been 
recalculated and the sixth degree polynomial perturbative corrections to its absolute minimum have been determined. 
Some of our conclusions have been compared with previous results obtained by other authors. 

To overcome some computational difficulties in the full characterization of the stratification of the orbit space, 
difficulties originating from the high dimensionality of the problem, we had to devise a sophisticated procedure, 
allowing to obtain rational parametric equations also for the higher dimensional singular strata. The procedure is 
absolutely general and is reminiscent, but in a different context and in higher dimensions, of classical techniques for 
parameterizing plane algebraic curves |23] . 



^The inclusion of terms of 6-th degree, depending only on pi,P2 and ps would be useless for splitting the degenerate minima, 
so these terms will be neglected, with no loss of generality. 
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APPENDIX A: CALCULATION OF MOLIEN FUNCTIONS 

In this Appendix, we shall explicitly calculate the two integrals, /i and I2, appearing in the expression of Mq^t]) 
(see (§)): 

1 Z"^" . , f^"" , r . ^,^ f^"" da 



1 / d02 / sinOde 



327r3 7o 7o 7o J„ /i(cosx,a)' 

I2 = 



1 Z"^" ., /■^'^ ., Z"' sin0d0 



(Al) 



where 



167r2 7o 7o 7o /2(cosx)' 

2 
/i(cosx,a) = n [(l-^e'^''^+"^)(l-^e'^'='^-«))" 

fc=-2 
2 

/2(cosx) - n [{^-Ve''''){l + ve'''')] 



(A2) 



fc=-2 



and cosx is defined in (|10[). 

Using as new integration variables 02, u = cos((/)i + (^2) and v = cosx, the integral over (f>2 can be immediately 
performed and one remains with 



I f ^ , f dv f du 



/i = 7—^ / da 



27r2 7o J-i fi{v,a) Jv [1 + u)Vl 



u^ 



= ^ r da t ^^ ^^^ (A3) 

= ^ r da r d^.'"-'^^"^ 



and, analogously. 



After defining 



Stt^ 7o Jo /i(cosx,a) 



1 f^- sin^(x/2) 
T^ Jo /2(cosx) 



(A4) 



Cl = e^^ C2=e'^", (A5) 



the remaining integrals can be conveniently thought of as integrals on unit circles of complex planes of the variables 
Ci and C2 

h = -^^Au Wr 7t/'^Ci^i(Ci,C2,^) 

8 7r^?7 / (C2 ~V)[C2V~ 1) / 

(A6) 



^^- 477^4^31) /^CrF2(Cr,.), 



where 
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Ct (Ci - 1) 



(C? - f ) (C? - C. V) (Ci - f ) (Ci - C. .) (Ci - f ) (Ci - c^) (c? - 1) (c? - ^) ^^^^ 

p r^ ^ C^ (Ci - 1)' 



(C?-^2)(Cf-ry2)(c?-J,)(c4-J,)' 

Using the Cauchy theorem and recaUing that \ri\ < 1, the calculation reduces to a computation of residues at the 
internal poles to the unit circles. So, one ends up with the following expressions for /i and I2' 

1 + 3 77** + 2 77!" + 3 77I2 + ,;20 



2 (1 - 772)3 (1 - 774)2 (1 - 778) (1 - 7; + 772)' (1 + 77 + 772)' ' 



(A8) 
1 — rj'^ + rf 



2(1- 7/2)3 ^^ _ ^4^ (1 - 7/ + 7/2) (1 + 7/ + 7/2) ' 



jk [-36 a' 6' + 144 6' 6 6' - 180 a' 6' 6^ + 72 6' ^3^ + 36 6' 6'- 
108 Ci' 6 6' + 72 6' 6' -9 a' 6' + 18 6 6' + 80 a' 6' a -144 666' C4- 

24 6 6' 6 - 48 6' 6' - 44 6' 6' 6 + 144 6 6 6' 6 - 72 6 6' 6+ 
1056 6' 6 6 + 144 6' 6' + 96 6' 6' 6 + 24 6' 6 - 384 66'+ 
24 01 6 (-86' 6 + 126 6 + 6' -166) -96 02 6 (66-46) + 

16 03 (6' + 246) -384 04]; 



(A9) 



h = 

The sum of h and I2 gives to the expression of Mc{ri) reported in (|1 1|) . 

As a byproduct of the calculation, one finds also the following expression for the Molien function of the subgroup 
SO3 of G: 

MsoM^^J dcf,, J d4>2 j d0sin0/i(cosx,O) 

= - A ■ 4!i 2T /dCiFi (6,1,77) 

4 7r7 77^ (1 — 7/^) J 

__ 1 + 774 + 7/8 

= (1-772)3(1-^3)4- 

APPENDIX B: SYZYGIES 
1. Explicit form of the syzygies for the group G 

The ten elements of the non minimal IB {6 0} defined in (|l4| , [l5| ) are related by the following syzygies: 

01^ = Tfe[36 6'6-i446'66 + i80 6'6'6-726'6-366'6' + i086'66'- 
72 6'6' + 9 6'6'- 1866' -806^66 + 1446666 + 2466'6 + 48 66'+ 
144 6' 6 6 - 144 6 6 6 6 + 72 6 6' 6 - 288 6 6 6 - 144 6 6' - 96 6' 6 6- 



(Bl) 



24 6' 6 + 192 6' + 24 01 (8 6' - 12 6 - 6) 6 + 96 02 (6 6 - 4 6) - 16 03 6] ; 

diO^ = -i(666 + ^i66-2026 + 2036); (B2) 



(B3) 
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72 C3 



yig [-36^1" 6' + 144 Ci' 6 6' -180 a' 6' 6' + 72^1 6' 6' + 36^1' 6'- 

108 Ci' 6 6' + 72 6 6' 6' - 9 Ci' 6' + 18 6 6 6' + 144 Ci' 6 ^4 - 576 ^i" 6 ^3 ^4+ 

720 Ci^ 6' 6 ^4 - 288 6' 6 a - 64 ^i" 6' ^4 + 288 a' 6 6' ^4 - 288 6' 6' ^4+ 

12 6' 6' ^4- 72 6 6' a -320 a' 6^4' + 576^1 6 6 a' +48 6 6' a' + 192 6^4' + 

144ei'6e5-576a4e2e3e5 + 720ei'6'^3a-288e2'e3C5-288ei"e3'e5 + 

576 Ci' 6 6' ^5 - 288 6' 6' 6 - 36 ^i" 6' Cs - 72 ^2 6' Cs + 256 ^i' 6 ^4 ^5 + 

1440 6 6' ^4 ^5 - 960 ^3 U^ Cs + 576 ^i' 6 ^5' - 576 ^i 6 6 Cs' + 432 ^i 6' a'- 

4800 a ^4 C5' - 576 ^3 Cs' + 96 ^^ ^32 ^^ + 24 fi 6' Ce + 384 ^^ ^3 ^^ ^^ _ 1152 ^ ^ ^ a_ (B4) 

192 6' ^4 Ce - 384 ^^ ^3 ^^ ^^ _ 96 ^2 ^^ ^^ _ 334 ^^ ^^ ^^2 ^ ^^gg ^^ ^^2 ^ ^^gg ^^ ^^2_ 

24 01 (8 6' 6' - 12 a 6 6' - 6 6' - 32 a' a a + 48 a a Cs + 4 a' a + le a a ^6+ 
64 a ^6 -64 a ^6) + 



(36 a' - 144 a^ a + iso a' a - 72 a' - 36 a' a + 108 a 6 a - 72 a' a- 

+ 72aac4 

2 000 t 2 , 



3 a' a' - 18 a a' - 80 a' a + 144 a a a + 72 a a u + 48 a' + 144 a' a 



144 a a ^5 + 168 a a Cs - 48o a a - 336 a' - 288 a^ a + 288 a a) + 
16 03 (a a^ - 4 a a + 96 a' a - 144 a es - 16 a 6 + 24 a a) - 384 e^ a] ; 

2^6 [- 144 a' 6' + 144 a' a a' + 1008 a' a' Cs' - 1872 a' a' a' + 864 a' a'+ 

108 a' a' + 144 a' a 6' - 1188 a' a' 6' + 936 a' a' - 144 a' 6 a' + 288 a' a^- 

9 a' 6' + 18 a 6' + 1728 a' a ^4 - 6912 a' a a a + 864o a' a' 6 a- 

3456 a a' 6 a - 1408 a' a' u + 5568 a' a 6' ^4 - 5184 a a' 6' ^4 + 416 a' a' a- 

1296 a a 6^^4-24 a a^e4- 3840 a^ a ^4' + 6912 a' a 6^4' + 960 a' a' a'- 

576 a a' ^4' - 48 a' a' + 2304 a a a' + 1728 a' a a - 6912 a' a 6 a+ 

8640 a' a' 6 a - 3456a a' 6 a - 2304a' a' Cs + 4032a' a 6' a- 

I728aa'e3'a-i584aac3'a-72aa^e5 + 3072a'ac4a + i4976a'a'e4a+ 

3456 aa'aa + 288 a'aa- 11520 aaa'a + 6912 a'aa'- 6912 a'ae3a'+ 

4032 a' a' a' + 1728 a a' a' + 144 a' a' - 57600 a a ^4 a' - 6912 a a ^5'+ 

864 a' a Ce - 3456 a' a 6 a + 4320 a' a' a a - 1728 a' a a - 480 a^ a' ^6+ 

3744 a' a a' a - 1728 a' a' a + 408 a' a' Ce - 144 a a' a + 24 a^ a+ 

2688 a' a a a - 10368 a a 6 a ^e - 1344 a a' a a - 384 a a' a - 1152 a' a es a- 

3456 a a 6 a Ce + 576 a a' Cs a - 14592 a a ^5 a - 3456 a a' ^e - 5376 a' a a'- (b5) 

2304 a a Ce' - 768 a' a' + 18432 a a Ce' + 18432 a a ^6^ + 6144 a'+ 

8 01 (36 a' a- 144 a^ a 6 + 180 a' a' a -72 a' a -132 a' 6' -36 a' a 6'+ 

240 a a' + 1296 a' a a - 1872 a a 6 ^5 - 72 a a' a - 864 a a ^5 - 144 a a'- 
96 a' a a - 864 a a a - 96 a' a - 2304 a a a + 2304 a a a + 768 a') + 

96 02 (36a'-i44a'a + i8oa'6'-72aa'-36a'a + io8a'6a-72aa'a+ 
5 a' 6' - 30 a a 6' - a a' - 80 a' a + 144 a' 6 ^4 + 72 a' a a + 48 a a'+ 
144 a^ a - 144 a' a Cs + 136 a' a Cs + 48 a a Cs + 4 a' a - 48o a a Cs- 
336 a a' - 288 a' a + 288 a a a + le a 6 a + 64 a a - 64 a a) + 

16 03 (36 a' - 144 a^ a + isoa'a' - 72 a' - 36 a' a + io8a'a6 - 72a'a+ 
13 a' 6' - 6 a a' + es' - 80 a' a + 144 a a a - 24 a a a + 48 a'+ 
1296 a' a - 1872 a a es - 72 a a a - 480 a a - 336 a' + 192 a' a - 48 a a) - 

512 04 (3a' + 9a +6)]; 



360 a' a' - 3 a' a' + 54 a a' + 3 a' - 80 a' a ^4 + 144 a a a ^4 - 24 a a' a 



71 173- 



(B6) 



2. Explicit form of the syzygy for the the rotation subgroup of G 

The invariants Ai, . . . , As of the subgroup SO3 of G, defined in ([49|), satisfy the foUowing syzygy: 
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80 Ai^ \i^ + 24 Ai Aa^ M^ - 16 Az^ A4^ + 24 Ai Ag^ A42 + 48 A2 Ag^ A4^ + 48 A4''- 

96 Aa^ A3 A4 A5 + 32 Ag^ A4 A5 - 80 Ai^ Ag^ + 24 Ai Aa^ As^ + 16 Az^ As^ + 24 Ai Ag^ Ag^- 

48 A2 Aa^ As^ + 96 A4^ Ag^ + 48 A5'' + 96 Ai^ A2 A4 Ag + 96 Ai Az^ A4 Ag - 48 Aa^ A4 Ag- 

96 Ai Aa^ A4 Ag - 48 A2 Ag^ A4 Ag + 96 Ai^ A3 A5 Ag + 192 Ai A2 A3 A5 Ag - 48 Aa^ A3 A5 Ag- 

48 A3^ As Ag + 144 Ai^ Ag^ - 288 Ai^ A2 Ag^ + 72 Ai Aa^ Ag^ + 72 Ai A3^ Ag^-- 

288 A4^ Ag^ - 288 As^ Ag^ + 384 A4 Ag^ - 144 Ag'' - 96 Ai^ A3 A4 A7 + 192 Ai A2 A3 A4 A7+ 

48 Aa^ A3 A4 A7 + 48 A3^ A4 A7 + 96 Ai^ A2 A5 A7 - 96 Ai A2^ A5 A7 - 48 A2^ A5 A7+ (B7) 

96 Ai A3^ As A7 - 48 A2 Ag^ As A7 - 576 Ai^ A3 Ag A7 + 1152 As Ag^ A7 + 144 Ai^ A7^+ 

288 Ai^ A2 A7^ + 72 Ai A2^ h^ + 72 Ai A3^ A7^ - 288 A4^ A7^ - 288 As^ A7^- 

1152 A4 Ag A7^ - 288 Ag^ A7^ - 384 A5 A7^ - 144 A7^ - 144 Ai'' Ag + 108 Ai^ A2^ Ag- 

18 A2'* As + 108 Ai^ A3^ As - 36 A2^ A3^ As - 18 A3'' Ag + 144 Ai A4^ As + 

144 Ai Ag^ As - 288 A2 A4 Ag As - 288 A3 As Ag As - 144 Ai Ag^ As + 288 A2 Ag^ As + 

288 A3 A4 A7 As - 288 A2 A5 A7 Ag + 576 A3 Ag A7 Ag - 144 Ai A7^ As - 288 A2 A7^ As + 

180 Ai^ Ag2 - 72 A2^ Ag2 - 72 A3^ Ag^) . 



APPENDIX C: EXPLICIT FORM OF THE P-MATRIX FOR THE GROUP G 



For 



(di,...,d9)==(2, 4, 4, 6, 6, 8, 8, 10, 12), (CI) 

the explicit form of the P-matrix related to the integrity basis {p} of the group G is the following 
Pi,i^2dip„ 2 = 1,..., 9, 
^2,2 = -3 (3pi3-9piP2-2p4 + 2p5), 

-^2,3 = 8piP3, 
P2,4=4 (2piP4+P6), 

4 
P2S, = o (14piP5 - 3pg - 2^7) , 

2 
P2fi = -^ (3p3P4 - 3p3Ps - 26pipg + 2ps) , 

4 
^2,7 = -o (2P3P5 - 17piP7 + 3pg) , 

2 

■P2,8 = -g (2P3P6-3P3P7 -32pips + 3p9), 

P23 = 4 (]53P8 + 5piP9), 

^3,3 = 16piP3, 

P^A = 24p6, 

-P3,5 = 8 (pg + 2p7) , 

-P3,6 = 4 (p3P4 + 3p3P5 + 2pipg + 2ps), 

^3.7 = 8 (2p3P5 +PlP7+Pii), 

P3.S = 4 (2p3pg + 3p3P7 + 4piPs+P9), 

^3,9 = 24 {psps +pipg), 

P44 = 3 (6pi^ - 6p2 -P3) P4, 

P4,5 = 2 (2pipg+p8), 

Pa,6 = -^ (4piP3P4 + 24p4ps + 18pi^pg - lSp2P6 - ipsPe + 2^9) , 

P4.7 = ^ (-36^1^ + 144pi4p2 - 180pi2p2' + 72 P2' + 36^1^3 - 108pi2p2P3 + 72p2'p3- 
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9pi^P3^ + lSp2P3^ + SOpi^Pi - U4piP2P4 - 24piP3P4 - 48p4^ - lUpi^p5 + 

U4pip2P5 - 72pip3P5 + 288p4P5 + 144^5^ - 96pi^pe + 288p2Pe + ll^p^pe + 
288pi2p7 - 288p2P7 + 32pip8 + I6P9) , 

Pi,S = 2 {'^PS^Pi + 48p4P7 + ISpi'^ps - 18p2P8 -3p3P8 + 4piP9) , 
■P4,9 = 3 {l2piPs + 6pi'^pg-6p2P9 -P3P9) , 
Pb.b = o (2P1^P5 +GP2P5 + 5p3P5 + 8P1P7) , 

^5,6 = ^ (-36pi6 + 144pi4p2 - 180pi2p2' + 72 P2' + 36^1^3 - l08pi^p2P3 + 72p2^p3^ 

Qpi^P3^ + 1&P2P3^ + 80pi^P4 - lUpiP2P4 - 24piP3P4 ^ 48 P4^ - 144^1^135 + 

lUpiP2P5 + 24piP3P5 + 480p4P5 + 144 P5^ - SOpi^pe + 336p2Pe + 88p3Pe + 
288 Pi^pr- 288 P2 P7 + 48 ^3 P7 - 32 pi pg + 16 pg) , 
1 



P5,7 = -^ (-36pi^ + 144pi''p2 - I80pi^p2^ + 72p2^ + 36pi^p3 - l08pi'^p2P3- 

2^2 I 10^ n^2 I Qf\ ^ 3 , 



96 

72p2^P3 - 9pi^P3^ + l8p2P3^ + 80pi^P4 - 144piP2P4 - 24piP3P4 - 
3 



^5,8 
-f5,9 
P&3 



48 p4^ - UApi^p5 + lUpip2P5 + 184pip3P5 + 288p4P5 + 912^5^ - 96pi2p6 + 
288p2P6 + 48p3P6 + 352pi^p7 - 96p2P7 + 160p3P7- 96piP8 + 16^9) , 

- (6p3^p5 +8piP3Pe + 96p5P6 + ^2pip3pr + A8p4pr + 2pi'^ps + 6P2P8 + 5p3Ps) , 

2 {-l2p3P4P5 + P3'^ pe + 24:P(i'^ + 2piP3Ps + 6p4P8) , 

1 

12 



{2 pi"^ P3 P4 + 6 P2 P3 Pa + 5p3^P4 + 54pi^p3P5 - Mp2P3Pb - 9p3^P5 +48piP4P5- 



Pf 



6,7 



6,8 



24piP3P6 + 48p4P6 + 144p5P6 + 96p4P7 + 12p3P8 + 8P1P9), 

4 (-36pi" + 144pi5p2 - 180pi3p2' + 72pip2' + 36pi5p3 - IQ^pi"" P2P3+ 
96 ^ 

72piP2^P3 - 9pi^P3^ + l?>piP2P3^ + 80pi'^Pi - lUpi^p2P4 - 2Api^p3P4 - A8 pi Pi^ 

UApi*p5 + lUpi'^p2P5 - 72pi^p3P5 + 96p3^P5 + 288pip4P5 + 144pip52 - 96pi^p6 
288pip2P6 + l76piP3P6 + 1152p5P6 + 288pi3p7 - 288pip2P7 + I92pip3pr + 
576p4pj + 576p5pj - GApi'^ps + 96p2P8 + 80p3Ps + IGpipg) , 
1 
12 



(8piP3^P4 + 48p3P4P5 + 12p3^P6 + 96p6^ + 54pi^p3P7 - Mp2P3P7 - Qps'^P?- 



96piP4P7 + 24:piP3Ps + 72p4Ps + 216p5Ps + 2pi' p^ + 6p2P9 + 5P3P9) , 
P6,9 = -^ (2p3^P4 + 2Ap3P4Pe + l8pi'^P3Ps - 18p2P3P8 -3]53^P8 + 2ApiP4Ps + 4piP3Pg + 

12p4P9 + 6OP5P9), 
^7,7 = ^ {2pi'^p3P5 + 6p2P3P5 + 5p3'^P5 + ^2piP5'^ + 24:P5Pe + 8piP3P7 + 48p5P7) , 

■P7,8 = 7 (l2piP3^P5 + 36p3P4P5 " 36p3P5^ + 2pi^p3P6 + ^P2P3P& + 5p3^P6 +48piP5P6 + 

6 

24 P6^ + 6p3^P7 + 72^7^ + 8piP3P8 + I20p5ps) , 
P7,9 = -12piP3P4P5 + 4piP3^p6 + 12^3 P4P6 ^ 12p3P5P6 + 2ApiPe'^ + 2p3^ps + 24prps + 
24p5P9, 

-P8,8 = -pj (2pi^P3^P4 + &P2P3^ Pi + 5p3^p4 + 54pi^p3^p5 - 54p2P3^P5 - 9p3^P5 + 

192piP3P4P5 + 24piP3^P6 + 96p3P4P6 " 288p3P5P6 + 12^3^ P8 + 576p7P8 + 
8piP3P9 + 192p5P9), 

•^8,9 = -^ (4piP3^P4 + 12p3^P4^ + 18pi^P3^P6 - 18p2P3^P6 - 'ip3^pe + A8piP3P4P6- 
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R 



9,9 



84p3P6^ -2Ap3PiPs + 84ps^ + 2p3^P9+ A8 pe pg + SA pj pg) , 

3 (Gpi^Ps^Pi - GP2P3^P4: -~P3^P4 + 12piP3^P4^ + 288pi P3^ P4P5 " 1152^3^4^5^- 
24p3^P4P6 - 288piP3P6^ + 1152p5P6^ + 144pi^p3P4P7 - 432p2P3P4P7 - T^P'i^PiP'j 

1152p4P7' + 108pi%8 - 432pi4p2P8 + 540pi2p2'p8 " 216p2'p8 - lOSpi^paPs + 

324pi^P2P3P8 - 216p2^P3P8 + 27pi^P3^P8 - 54p2P3^PS - 240pi^P4P8 + 

432 Pi p2 P4 P8 + 72 Pi p3 Pi p8 + 144 pi' ps + 432 pi^ p^ ps - 432 pi p2 pz Ps + 
2l6pip3P5ps - SGApiPzps - 432p5^p8 - 86Api'^pjps + 864p2P7P8 + U'ipi'^ P5P9 - 
4:32p2P5P9 - 72p3P5pg + 288piprpg) ■ 



APPENDIX D: EXPLICIT FORM OF THE GENERATORS OF THE ISOTROPY SUBGROUPS OF SO3 

In this Appendix arc collected the definitions of the groups and of the elements and the generators of the subgroups 
of SO3, involved in the body of the paper, in the tables and in the figures. 
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a. Generators or group elements 
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C2. = ( -10 

1 



G 



ix 




R.(0) = 



cos(0) sin((/)) 

— sin((/i) cos((/)) 

1 




1 

R,(0) = ( cos(</)) sin((/.) 
— sin(0) cos{(j)) 



(Dl) 



h. Groups 



By SO3 and O3 we denote the proper and, respectively, the complete 3D rotation group and, by Oj, a 



x,z the following O2 subgroups of SO3: 

{C2X R-z(0)}o<0<27r- 



Of = {Rxi(f>)}0<^<2^ U {C2.Rxi(f>)} 



0<4><2-!7 



and Of = {Rz{<P)}o<4><27r U 
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TABLE L Relations defining ID geometrical strata, S^'^"^^ = Vl^(i"^) in orbit space. For 2 < i < 9, p; = p^/{p'f'^^), di 
denotes the degree of the polynomial Pi{x) in the order parameter x and the values, p^^' ' , taken on by p on W^^' ' are listed 
in the columns from 2 to 9. 
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P5 
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1/2 
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1/6 


1/6 


1/6 


1/6 


1/6 
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1/2 
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1/3 





1/3 

















5 


1/2 
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TABLE II. Parametric equations defining tlie geometrical strata, W^^' ' (A is listed in the head row) in orbit space. For 
2 < i < 9, Pi = pi/{p^^' ) and di denotes the degree of the polynomial pi{x), in the order parameter x, and e — ±1. 

p\A U 2 3, 4 5 



0/12 



P2 




(2 + 0)/6 


P3 







P4 


(2 


-C)'(i + 


P5 







P6 







Pr 







P8 







P9 







^ range 




]0,^[ 



1/2 


{2 + e)/6 


1/2 


e 


? 


e 


1 


2-2C 





(2-0'(l+C)/12 


^ 








e'(i+c)/i2 


^ 








(2-C)(l + 0CVl2 


^ 








(i+c)eVi2 


^ 








(2-C)(l + 0CVl2 


^ 








(2-0'(l + C)CVl2 


^ 





o,i[ 


]o,.[ 


]o,i[ 


]il[ 



TABLE III. Possible symmetry strata, S^'''^' for D-wave driven pairing in isotropic space (d denotes the dimension of the 
stratum and r is an enumeration index). From left to right, the columns refer to the phase reference numbers according to 
our ((d, r)) and ref.[12] (A/") classification, the order {\H\) and the the residual symmetry group, H, the complex coordinates 
{zj — Xj + ix^+j, j = 1, ... ,5) of the H-invariant superconducting vector order parameter and corresponding values of the 
partial wave amplitudes (-D2, • • . , D-2). In columns five and six, the i's are real parameters, while the w's are complex ones. 
The definitions of the group elements are recalled in Appendix D. 



(d,r) 


N 


\H\ 




H 


(2:1,..., Z5) 


V2{D2,Di,Do,D-i,D-2) 


(1.1) 


II 


oc 


(C2 


^r)x{R,(0)Ui(20)}^ 


(it, -t, 0,0,0) 


(2t, 0,0, 0,0) 


(1,2) 


VIII 


00 




o|®(r) 


(0,0,0,0,t) 


(0,0,tv^,0,0) 


(1,3) 


IX, X 


16 




(C2^.r,C4.Ui(7r)> 


(D,t, 0,0,0) 


(-t,0,0,0,-t) 


(1,4) 


XI 
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(C2 


,,C2ar,C3iUi(4,r/3)) 


(0,-«,0,0,t) 


(«,0,t^,0,«) 


(1,5) 


I 


00 


(C2 


,r>x{R,(«Ui(-0)}^ 


(0,0,t,-«,0) 


(0,0,0,2 it, 0) 


(2,1) 


V 


6 


(C2a:-r,C3,Ui(47r/3)) 


(itl,-ti,t2,-it2,0) 


(2ti,0,0,2it2,0) 


(2,2) 


IV 


4 




<C2,r,C2,Ui(7r)) 


(0,0,ti ,1*2,0) 


(0,iti+it2,0,iti-it2) 


(2,3) 
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8 




(C2=.,C4.r> 


(0,itl,0,0,t2) 


(-iti,0,t2V2,0,-iti) 


(2,4) 


VII 
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(C2^.C2.,r> 


(0,ti, 0,0,(2) 


(-ti,0,t2V^,0,-ti) 


(2,5) 


VI 
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(C2xr,C4.Ui(7r)) 


(itl,t2, 0,0,0) 


(ti-t2,0,0,0,-ti-t2) 
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(C2,,C2,r> 


(«l,t2,0,0,t3) 


(ti-t2,0,t3V^,0,-ti-t2) 


(3,2) 
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(C2=.,C2.) 


(0,ui,0,0,-i;2) 


(-vi,0,V2V2,0,-vi) 


(4,1) 
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(CaxT) 


(«l,*2,t3,«4,*5) 


(tl-t2,it3 + it4,t5^,it3-it4,-ti-t2) 


(4,2) 


— 
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(C2.) 


(^1,^2,0,0,^3) 


( — iv-i—V2,0,V2^,0,ivi —V2) 


(6,1) 


-- 
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(vi,V2,V3,Vi,v^) 


( — ivi —V2,iv^+V4,v^^/2,iv^—V4,ivx —V2) 
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TABLE IV. Strata of dimensions > 3, bordering each stratum S^''''^' {d denotes the dimension of the stratum and r is an 
enumeration index). From left to right, the columns refer to the enumeration numbers, (d, r), of the lower bordering stratum, 
S^'^''^' , its orbit-type, [G''*''^'], the maximal conjugacy classes of subgroups of [G''*'''], [H]ai, the order, \Hm\, of Hm and the 
indices, {d',r'), of the corresponding strata, if any. The definitions of the group elements are recalled in Appendix D. 
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[G( 



d,r)l 



[Hh 



\Hm\ 



{d',r') 



(2,1) 


[{Ga.T, G3.(7i(47r/3))] 


[(G3.t/i(4V3))] 
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[{C2.T)] 
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(2,2) 


[{G2.T,G2.t/l(7r))] 


[{C2.T)] 


2 






[(C2.Ui(7r))] 
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(2,3) 


[{G2.,G4xT)] 


[(G2.,G2.T)] 
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[(G2.,G2z) 
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[(C4.T)] 
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(2,4) 


[{G2.,G2.,T)] 


[(G2.,G2.)] 
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[(C2.,G2.T)] 
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(2,5) 


[{G2.T,G4.C/i(^))] 


[(C2.,G2.T)] 


4 






KG4.C/i(^))] 
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(3,1) 


[{G2.,G2.T)] 


[{C2.T)] 


2 






[{C2.)] 
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(3,2) 


KG2.,G2.)] 


[{C2.}] 
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(4,1) 


[{C2.T}] 


[1] 
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(4,2) 


[{C2.}] 


[1] 
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4,1 
4,1 

3,1 
3,2 

3,2 
3,1 
3,1 

4,1 
4,2 
4,2 
6,1 
6,1 



TABLE V. Absolute minimum, ^^^l^ = — ai/(25), of a general, bounded below, G-invariant 4-th degree polynomial. 



$W(a,p) =«oP?/2 + E,=i 



«jPj, "1 



of the strata is the same as in Table fl% 



< 0, and hosting strata, S^''''^' , as functions of the coefficients a. The denomination 



Q range 



7^l : Max(0, -3q:o/2, -6Q3) < Q2 

7?.2 : —603 > a2 > Max(— Qo — 203,203) 

7^3 : -ao/2 <a2 < Min(0,2a3) 

7?.i3 : = 02 < Min(ao,a3) 

7^l2 ; Max(-3ao/2,0) < 02 = -603 

7?.23 : — ao/2 < Q2 = 2q3 < 

7?.i23 : Q2 = Q3 = < ao 



(d,r) 



ao + 2a2/3 

ao + a2 + 2q3 

ao + 2a2 

ao 

ao — 4a3 

ao -I- 4a3 

ao 



(1,4) 
(1,2), (1,3), (2,4) 

(1,1) 

(1,1), (1,4), (1,5), (2,1), (3,1), (4,1) 

(1,2), (1,3), (1,4), (2,3), (2,4), (3,2) 

(1,1), (1,2), (1,3), (2,4), (2,5) 

all, except (0,1) 
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FIG. 1. Possible phase transitions between bordering strata, connected, in the figure, by continuous sequences of one or more 
arrows. The group elements are defined in Appendix D. 




S(i.s) 
[{C2:,T) X {R.(0)Ui(-0)}^ 



5(1,4) 
[(C2^:C35Ul(4f),C2ar)] 



[(C2, 



S(l,3) 

,C4. Ui(7r),r)] 



5(1,2) 

[Of X (T) ] 



5(1.1) 
[(C2=.r) X {R,(0)Ui(20)}^ 




{e} 



FIG. 2. Localization of the absolute minimum of a fourth degree G-invariant polynomial, $'''' (a, p) = aop\/2 + ^ ,_ ajPj , 
Qfi < 0, as a function of its coefficients. For values of (0:0,02, as) in Tii or in 7?.2 or in TL^, the absolute minimum lies, 
respectively, in the strata S'^'*' or {S'^'^ , 5''^''^ , 5'^'^'} (degenerate minimum) or S''^'^'. For particular values of the a's, see 
Table livl 
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a2 



TZ2: ao > -2 0:2/3 



TZi : ao > —2 0:2 



"2=0 



-|> 



03 



7?.3 : ao > — a2 — 2 03 
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